arXiv:1502.05415v3 [nucl-th] 2 Jun 2015 


UT-15-05 
February, 2015 


Diquark bound states with a completely crossed ladder 

truncation 


Ryusuke Jinno, Teppei Kitahara, and Go Mishima 


Department of Physics, Faculty of Science, University of Tokyo, 
Bunkyo-ku, Tokyo 113-0033, Japan 


Abstract 

The Bethe-Salpeter equation in the diquark channel is investigated by employing the 
Dyson-Schwinger method together with the Munczek-Nemirovsky model. The novelty of 
our study is a resummation of completely crossed ladder diagrams in the Bethe-Salpeter 
kernel. These diagrams are enhanced due to their color factors in the diquark channel, 
but not in the meson channel. In our analysis, diquark bound-state solutions exist in the 
Bethe-Salpeter equation. 
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1 Introduction 

One of the greatest challenges of particle physics is to understand the dynamics and mass 
spectra of hadrons in terms of Quantum Chromodynamics (QCD), and a concept of diquarks 
- quark-quark bound states - has been playing an important role towards these goals. For 
example, diquark-quark models of baryon [1-3] are used to reproduce the mass spectra, the de¬ 
cay rates, the magnetic moments, the ratio of the axial-vector coupling to the vector coupling, 
the charge radius, and the processes of deep-inelastic scattering [4,5]. Another example is phe¬ 
nomenology of recently confirmed exotic hadrons A(3872), Zc(3900)^, Zc(4020)^, Z(4430)^, 
etc. [6-10]. It is proposed that a diquark-antidiquark description of such exotic hadrons can 
explain observed properties of those states, especially the narrow width and decay proper¬ 
ties [11-15]. 

Meanwhile, a functional method using the Bethe-Salpeter (BS) equation with the Dyson- 
Schwinger (DS) equation has been considered to have no diquark bound-state solutions. In 
the case of mesons, a consistent scheme which preserves the Nambu-Goldstone nature of pions 
is established for the extensively used rainbow-ladder (RL) truncation [16] and its extension 
[17-19], and the resulting BS eq. has bound-state solutions in both truncations. However, in 
the diquark channel of the BS equation, a certain kind of beyond-RL contributions have been 
found to remove the bound-state solutions [19-22], in contrast to the RL truncation [23,24]. 

In this paper, we investigate the differences among previous studies [19-22] which make 
the diquark bound-state solutions absent, and propose another approximation scheme in which 
diquark bound-state solutions exist. To this end we employ the same QCD model as Refs. [19, 
21,22], i.e., the Munczek-Nemirovsky model [25]. The simplicity of this model enables one to 
take higher loop diagrams into the BS kernel [19,21,22,26-32], leading the previous studies to 
find that a certain crossed ladder diagram gives a large repulsive contribution in the diquark 
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BS equation due to its enhanced color factor [19,21,22]. We point out that this enhancement 
is not peculiar to the diagram considered in these studies, but common among a class of 
diagrams - the completely crossed ladder diagrams. Then we perform a resummation of these 
diagrams adopting large-counting as a criteria for the truncation [33], finding bound-state 
solutions to the diquark BS equation. Although our calculation is based on a specific model 
with a specihc truncation scheme, we believe that our result is at least a counterexample of 
the previous arguments that the BS eq. has no solutions in the diquark channel. 

The organization of this paper is as follows. In Section 2, we introduce the DS eq. for 
quarks, and explain the model we adopt. Next in Section 3, we introduce the BS eq. for 
mesons and diquarks, and discuss the type of diagrams which give the leading contribution 
to the BS kernel. Our main analyses are presented in subsection 3.2 and 3.3. The numerical 
solutions of the BS eq. and interpretations of them are shown in Section 4, and the final 
section is devoted to conclusion. 

2 Dyson-Schwinger equation 

Dyson-Schwinger (DS) equations are functional relations [34, 35] which are derived from a 
functional change of variables in the path integral in a modern point of view. They include 
information on non-perturbative effects in principle, and provide us a useful way to study the 
infrared behavior of QCD. The DS eq. of the quark propagator takes the form 

S~^{p) = Z2{i:^ + mb)+T.{p), (1) 

S(p) = Ie^^S{q)t'^T,{q,p)D^,ik), (2) 

where S{p) is the dressed quark propagator, S(p) the quark self-energy, the dressed gluon 
propagator, T^ the dressed quark-gluon vertex, nib the bare mass of the current quark, the 
generator of SU(3)c, the gluon momentum k = q—p, g is SU(3)c gauge coupling constant, and 
Zi and Z 2 are the renormalization factors of the quark-gluon vertex and the quark propagator, 
respectively. The general form of the dressed quark propagator is given as 

Sip) = -^pA{p^) + Bip^). (3) 

The dependence on the renormalization scale C is omitted for simplicity. Here we use the Eu¬ 
clidian metric and the gamma matrices satisfy 7 j/} = 25 The dressed quark propagator 
is renormalized by the following condition: 

5’“Hp)lp2=c2 = i|> + mo(C), (4) 

where mo(C) is the current quark mass at the scale (. 

Equation (1) becomes a self-consistent equation of the quark propagator S(p) when F^ip, q) 
and Dfj,y{k) are modeled as some known functions of p, q, k and Sip). This method can realize 
large quark mass at low energy, which is considered to be a consequence of the dynamical 
chiral symmetry breaking [16,25,39,40]. 

In this paper, we adopt the Munczek-Nemirovsky model [25], 

g^D^yik) = (V - hK/k^) Gi27r)^6\k) 

= l5^yGi27r)^5\k), (5) 
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Figure 1: Solutions of the DS equation (1). The solid (dashed, dotted) line corresponds to the self¬ 
energy of (6) ((18), (26) with n = 6). The vertical axis is nondimensionalized by rj. The solutions for 
n = 3,4, 5 in (26) come between the red and yellow lines. See Appendix B for details. 

where G = rf /A is a dimensionful parameter which determines the scale of this model. In 
addition, the projection operator {S^y — k^ky/k'^) becomes due to the spherical symmetry 
of the model. In this model, the three-point vertex of gluon vanishes because of its derivative 
coupling property, but the four-point vertex still remains. Here and in what follows, we 
assume that the effects of the gluon self-coupling interactions are not crucial for the calculation 
of the bound-state mass. The validity of this assumption is not obvious especially when 
we consider large-A^ counting since vertex corrections due to these interactions are known 
to be enhanced in the large-counting [29,36-38]. We consider that the assumption is 
partially supported by the success for reproducing meson mass spectra in this setup as shown 
in Refs. [19, 21, 22, 25-30] and in Table 1 of this paper. We take the renormalization scale 
Q to be large enough (C S> 1 GeV), and we can set the renormalization constants to 1 and 
need not to specify the regularization scheme in the practical calculation since the model has 
no UV divergence. The dimensionful parameter r/ can be fitted to reproduce the observed 
mass spectrum well, and the value ry ~ 1 GeV is found to satisfy this requirement [25]. The 
5-function behavior in Eq. (5) is considered to be a realization of the property of hadrons 
dominated by the infrared dynamics of quarks and gluons inside them. Here the infrared 
region means the energy scale smaller than G ~ (0.5 GeV)^. 

In the Munczek-Nemirovsky model, the self-consistency equation becomes simple algebraic 
equations of A and H, and thus can easily be solved once the self-energy S(p) is specihed. One 
choice is the rainbow truncation. 


S(p) = 5(^)7^, (6) 

where the integration with respect to the loop momentum can easily be performed due to the 
5-function. The quark propagator functions A and B are evaluated by solving Eq. (1) with the 
self-energy Eq. (6), and the results are plotted in Eigure 1. We take the current quark mass 
to be mo = 0.012 x rj [GeV] in the plot, which is consistent with previous studies. With this 
choice, we obtain reliable masses of pion and rho mesons if we set ry ~ 1 GeV (see Table 1). 
More precise parameter fitting is possible but is beyond our scope; hence, we use this value 
for the current quark mass unless mentioned explicitly in the following calculations. 
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Figure 2: Diagrammatic expression of the homogenous Bethe-Salpeter equation in a diquark channel. 

3 Bethe-Salpeter equation 

Bethe-Salpeter (BS) equation has been the central tool to investigate bound states in terms 
of quantum field theory since its birth [41,42]. It is a functional relation between four-point 
correlation functions, and when we regard their residues as bound-state poles, the BS eq. 
becomes a homogeneous equation with respect to a certain function ~ BS amplitude. The 
homogeneous BS eq. takes the following form; 

nP,p) = 1-0^ K{p,p,k)s{k+)r{p,k)s{k_), (7) 

where p and P are the relative momentum and the total momentum of the system, respectively, 
T{P,p) is the BS amplitude and K{P,p,k), the BS kernel, is the qq scattering kernel which 
is two-particle-irreducible [43] with respect to qq lines, and we define = I ± P/2. The 
flavor, Dirac, and color indices are omitted in Eq. (7) for simplicity. Note that the way to 
contract these indices depends on whether the bound state is meson or diquark. Figure 2 is a 
diagrammatic expression of the homogenous Bethe-Salpeter equation in the diquark channel. 

In the Munczek-Nemirovsky model, the BS amplitude is a function of only P because the 
momentum transfer by gluons is zero, and its general form for the pseudoscalar meson channel 
is 


r(P) = 75[/i(^')+if’/2(P')], (8) 

while for the vector meson channel 

4 bMP") + , (9) 

where (A = —1,0,-|-1) is the polarization four-vector and = [ 7 ^, 7 j^]/ 2 . Substituting 
these into Eq. (7), one obtains a 2 x 2 matrix eigenvalue equation 

and a similar equation for the vector channel. Note that there is no remaining integral due 
to the (5-function in Eq. (5). In order for Eq. (10) to have a nontrivial solution, the following 
condition: 


det [/ - H{Pb] = 0, (11) 

with I being the 2x2 unit matrix must be satisfied. Conversely, if Eq. (11) holds for a certain 
value of P^ = —m? < 0, then the BS eq. has a bound-state solution and the pole mass of 
the state is given by m. We treat Eq. (11) as an equation of one variable P^, and investigate 
whether it has a solution or not. 

To evaluate the LHS of Eq. (11), one has to specify the diagrams taken into the kernel 
K. In the case of mesons, it is known that a consistent scheme for the DS and BS kernels 
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Figure 3: A diagrammatic expression of the Ward-Takahashi identity preserving combination of the 
quark Dyson-Schwinger eq. and the meson Bethe-Salpeter kernel in the rainbow-ladder approximation 
(0(G)) and next-to-leading order approximation (0(G^)). 



exists which automatically preserves the axial-vector Ward-Takahashi identity and keeps the 
Nambu-Goldstone nature of the resulting pion mass [16,17,19]. In this construction, the BS 
kernel is determined to be the sum of the terms with one of the quark propagators in the 
self-energy (2) replaced by 

s(p) ^ -s+r(p)s., (12) 

and with all the other quark propagators in the left (right) to the replaced one substituted 
by 5+ (S-). Here we defined S± = S(iLP/2). Following this scheme, one of the frequently 
adopted methods is to determine the self-energy (2) by the coupling expansion, and determine 
the diagrams in the BS kernel by Eq. (12). 

On the other hand, no such symmetry is known for the diquarks. This means that the BS 
kernel does not necessarily have a correspondence to the DS kernel like Eq. (12). 

In the subsection below, we examine the coupling expansion in both meson and diquark 
channels following previous studies. This enables us to infer the kind of diagrams - completely 
crossed diagrams - which have the largest color factors at fixed orders of G. Adopting a 
truncation scheme which takes these diagrams into account, we will see that the coupling 
expansion is unreliable in determining the existence of diquark bound-state solutions. Then 
we show that the resummation of all the completely crossed ladders, rather than the coupling 
expansion, leads to diquark bound-state solutions. 


3.1 Bound-state solutions up to O(G^) 

We begin with 0(G). The truncation of this order is called rainbow-ladder (RL) truncation. 
Applying the prescription (12) to the quark self-energy shown in Eq. (2), we have the following 
BS eq. for flavor non-singlet mesons 

r(p) = -Gj^s+r(P)s.j^. (13) 

The BS kernel of Eq. (13) is equal to the one obtained by the coupling expansion at this order. 
The diquark BS eq. with the leading contribution to the BS kernel is 

r^(P) = -lG7^5+r"(P)5_7;.. (14) 

Here some comments are in order for the diquark BS equation. The diquark BS amplitude 
has color indices due to its color non-singlet nature. For example, the diquark BS eq. with 
this truncation becomes 

r,,(p) = -lG(PU(nnf.s+rki(P)s^7^, (i5) 


5 















where T denotes transposition, i,j, k, I are color indices running from 1 to 3, and a is an adjoint 
index running from 1 to 8 . The Dirac matrices *S'+, 7 ^ and the color matrix are transposed 
in order to contract two quark lines correctly. The scalar diquark BS amplitude is expressed 
as 




/f (P2) + 


C = eijkr’‘iP)C, 


(16) 


where Cijk is the antisymmetric tensor and C = 7274 is the charge conjugation matrix. It 
satishes 


C2 = -l, ^7^6- = 7 m, CSlC = -S.. (17) 

Substituting Eq. (16) into Eq. (15) and using Eqs. (44) and (17), we obtain Eq. (14). 

We can set up the determinant equations for mesons and diquarks following the procedure 
( 8 )-(ll) using the quark DS eq. with the self-energy ( 6 ). At this order both BS eqs. have 
bound-sate solutions, and we summarize the resultant masses squared in the first column of 
Table 1 as vr, p and d{A). In addition, we plot the behavior of the LHS of Eq. (11) for diquarks 
in Figure 4, which actually shows that the line has an intersection with the horizontal axis. 
Note that the expression of the diquark BS eq. (14) has the same form as the meson one (13) 
except for the factor In general, the orders of 5+, 5_, 7 ^ appearing in the meson and 
diquark BS eq. coincide at any order. 

Next let us see 0{G‘^). The quark self-energy is given by 

S = G'jf.S-ffj, + (18) 

With the substitution (12), the meson BS eq. becomes 


T = - G^^S+TS.-i^ 

- |G^7^5+ [r5_7i.S'_7^ -k 7 i.S'+rS'_ 7 ^ -k 7j.5+7^S'+r] (19) 

Notice that there are three 0{G‘^) terms, since the replacement Eq. (12) is applied to the three 
quark propagators in Eq. (18). The diagrammatic expression of the consistent combinations 
of the self-energy part of DS eq. and the meson BS kernel is shown in Figure 3. Again, at this 
order, the BS kernel of Eq. (19) is found to be the same as the one obtained by the coupling 
expansion. The diquark BS eq. at this order becomes 


T^ = - iG7M5+r"5_7M 

- [r^S_7,5_7^ + 57,5+r^5_7^ + 7 , 5 + 7 ^ 5 +r* 


S-'yy. 


( 20 ) 


Solving the determinant equation (11) with the quark propagator calculated by Eq. (1) 
using the self-energy Eq. (18), one finds that the bound-state solutions for both pseudoscalar 
and vector mesons still exist. The resulting masses squared of these particles are listed in 
the second column of Table 1 as vr and p. In both channels, the difference between the RL 
truncation (13) and the 0{G‘^) truncation (19) is less than 10 %, and one can see that these 
masses are rather stable against the coupling expansion. In contrast, one finds that the bound- 
state solutions for diquarks disappear [19]. To see the disappearance of the solution, the value 

Since color-symmetric (6) diquarks do not form bound states even in this truncation, we consider only 
color-antisymmetric (3) diquarks in this paper. 
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Table 1: The values of bound-state mass squared, i.e., —P^ in rf [GeV^] unit obtained by solving the 
Bethe-Salpeter equation for pseudoscalar meson (tt), vector meson (p), and scalar diquark (d) channels. 
Hyphen means that there is no bound-state solution and we list all the solutions with — 6 < < 6 if 

the BS eq. has more than one solution. See the main text for the explanation of Scheme A, B in the 
diquark case. 



G 


G^ 


G5 

G^ 

TT 

(0.136)^ 

(0.140)^ 

(0.142)^ 

(0.142)^ 

(0.142)2 

(0.142)2 

p 

(0.724)^ 

(0.793)^ 

(0.824)^ 

(0.849)^ 

(0.859)2 

(0.842)2 

d{A) 

(1.14)^ 


(0.88)^ 

(1.97)^ 

(0.86)2 

-(0.58)2 

(1.09)2 

(1.16)2 

d{B) 

(1.14)^ 


(1.22)^ 

(0.31)^ 

(0.98)2 

-(0.53)2 

(0.34)2 

(1.18)2 

-(0.99)2 

(1.96)2 


of the determinant (11) calculated using Eq. (20) is plotted in Figure 4, where the plotted line 
has no intersection with the horizontal axis. 

Here note that there is ambiguity about which truncation scheme in the quark self-energy 
should be used for diquarks, Eq (6) or Eq. (18), because unlike the meson case, no symmetry 
relates the quark DS eq. and the diquark BS kernel. However, we confirmed that the disap¬ 
pearance of the diquark bound-state solution happens with both quark propagators; therefore, 
we do not discuss the quark DS eq. here. In fact, as we show later, all the results below are 
qualitatively independent of the choice of the truncation used in the quark DS eq. as far as 
we use the quark propagators listed in this paper. 

Where does the qualitative difference between mesons and diquarks come from? Recall 
that the bound-state solution is present for mesons at this order. The only difference in 
Eq. (19) and Eq. (20) is the factors on each term, especially that of the second term of O(G^). 
This term comes from the crossed-ladder diagram, which can be considered to give a repulsive 
effect. In order to check that, we solve the quark BS eq. taking only the crossed-ladder term 
in 0{G‘^) into account 

r" = - ^GS^5+7.5+r'=5-7M^-7., (21) 

and obtain qualitatively the same result. 

Note that the enhancement of the crossed ladder diagram is not peculiar to O(G^). To see 
this, let us see the color factors of (crossed) ladder diagrams at 0{G'^). The ladder diagram 
for the meson channel has a factor (see Appendix A) 






whereas the completely crossed ladder diagram has 

^ 0-2 ^ ^ ^ ^< 3-71 —1 ^< 3-1 ^^ 3-2 . . . -j^^n 


^n/2n. 

(22) 

• lV/2”, 

(23) 


in the large-limit. The enhancement of the ladder diagrams and the relative suppression 
of others including the completely crossed ladder diagrams provide one ground for the RL 
approximation. On the other hand, the ladder diagram for the diquark channel has 




To,n-l 


^Ta2^Tai 


e/2" 


(24) 
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Figure 4: The values of the det{I — H) in the diquark channel with Scheme A (left panel) and B (right 
panel). See the main text for the explanation of Scheme A and B in the diquark case. 

whereas the completely crossed ladder diagram has 

. . . ^a„_i^a„g^Tai^Ta2 . . . ^Ta„_i^Ta, ^ eN'^-^/2^, (25) 

where e is a symbolic expression of the antisymmetric tensor Cijk- 

What we find from these calculations is that the crossed-ladder contributions are not 
perturbations in the case of the diquark BS equation. Although the argument based on 
the large-A^ limit may not be directly applicable to the real-world QCD with N = 3, the 
enhancement found above indicates the possibility that an inclusion of higher order terms of 
this type changes the result again. 

3.2 Bound-state solutions beyond 0{G'^) 

First we investigate the effect of the completely crossed ladder diagrams on the meson BS eq. 
beyond O(G^) for comparison with that on the diquark BS equation. In order to realize both 
the inclusion of the completely crossed ladders and the symmetry preserving scheme explained 
in Eq. (12), we choose the self-energy of the quark DS eq. to be 

S = - Y^i-cr (26) 


with 



N - 


N 


2"+i 




j_N n-l' 

n) ’ 


(27) 


and Figure 5 is a diagrammatic expression of it. Then the replacement (12) determines the 
BS kernel. At 0(G”), there are 2n — 1 corresponding terms in the BS equation, and we can 
solve the BS eq. in exactly the same way as in the previous subsection. We put the actual 
calculation in Appendix B. The results up to 0{G^) are summarized in Table 1, which shows 
the effect of the completely crossed ladders is insignificant in the meson channel. 

Next we investigate the diquark BS eq. with the completely crossed ladder diagrams. It 
has a form 


= - j;(-G')" (f)" 4°)7m^+ • • ■ • • • 5_7;., 


n 


(28) 






























n 


n 


Figure 5: Diagrammatic expression of the quark self-energy of the completely crossed gluon correction 
(26). 

with 


4°) = ^ [(-l)yiV + 1)^ - ^ (iV - . (29) 

There are some options for the quark self-energy similar to what we have done in the 0{G'^) 
case. In order to investigate the truncation dependence, we try two schemes, Eqs. (6) and (26), 
for the quark self-energy, and call them Scheme A and B, respectively. Scheme A is chosen 
to include the leading-order contribution in the ’t Hooft limit [33], i.e. the large-A^ limit 
with G = G/N, which is consistent with the criteria used for the resummation of completely 
crossed ladder diagrams in the next subsection. Scheme B is the same one for mesons at 
each order of G. The determinant (11) for these diquark BS eqs. are plotted in Figure 4 
up to 0(G®), and the bound-state masses squared (if the solution exists) are summarized in 
Table 1. When the BS eq. has more than one solution with —6 < p2 

< 6, we list all of them 

in Table 1. As discussed in the previous subsection, the bound-state solution present in the 
0{G) truncation (the RL truncation) disappears at 0{G‘^). However, it revives at 0{G^) and 
somehow persists up to 0{G^). In Figure 4, one finds that the behavior of the determinant (11) 
at time-like region changes qualitatively as the order of G increases, especially that the sign of 
the gradient changes alternately. We try taking account of other diagrams in the diquark BS 
kernel by imitating the meson BS kernel constructed below Eq. (26), and obtain qualitatively 
the same result. 

From the enhancement of the color factors discussed in Eqs. (22)-(25) and the results 
displayed above, one possible origin of instability against the coupling expansion in the diquark 
BS eq. is the completely crossed ladder diagrams, which have the largest color factors at fixed 
orders of G. In the following, we provide a way to take all these diagrams into account. 

3.3 Resummation of the completely crossed ladder diagrams 

As we see in the previous subsection, the coupling expansion seems not to be valid in the 
diquark BS kernel due to their enhanced color factors in the completely crossed ladder dia¬ 
grams. Since the completely crossed ladders are found not to be perturbations, we perform 
a resummation of these diagrams in the diquark BS kernel. Note that all of these diagrams 
are the leading contribution in light of the’t Hooft limit [33], namely, the large-A limit with 
G = G/N. Therefore, our resummation can be regarded as a truncation with respect to the 
’t Hooft limit. 

The basis for the BS kernel is spanned by 16 elements, which we choose as follows: 

H = {6/11 = 1,2,..., 16}, (30) 
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where 


61 = 1 (g) 1, 

62 = 1 < 8 > 

63 = |^(g 1 , 

bi = f 


^5 = 7m 7m> 
fee = 7 m ® ^7m> 
fe? = I^7 m ® 7m> 
fes = -^7 m ® -^7mi 


feg = 7 m7i/ ® 7M7i^, 

felO = l^llu ® tliTiv, 
fell = fliTiv ® liiHv, 
fel2 = tliTiv ® fliilv, 


fel3 = liilvlp (g Iplvlp, 
fel4 = IpHvlp <g tlplvlp, 
fel5 = tlplvlp g 7 m7j^7m: 
feie = tlplvlp g f’lplulp 


(31) 


Other elements than those shown here (e.g. 7 /x 7 i/ 7 p 7 o- g Iplulpla) can be expressed by the 
ones shown here, as proven in Appendix C. Here bj = uj ® Ij is understood as (bj) 0 ^ 02,^1132 — 
{ui)oia 2 {h)i 3 ii 32 i where ai, a 2 , /3i, /32 are Dirac indices. Note that the Lorentz indices are 
contracted between the upper {ui) and lower (//) halves of the basis, since both are connected 
by gluon lines, which are effectively proportional to Also note that the only momentum 
appearing above is P, the total momentum of the diquark, since the gluon momentum vanishes 
in the present model. 

We dehne Kn as the completely crossed diagram in the diquark channel, with n gluon lines 
propagating between two quark propagators. Then our choice of BS kernel is K = X)^i 
(see Figure 6 ). Let us factorize Kn as Kn = {—G)^{j)^d^^Kn, where is the color factor 
dehned in Appendix A. The vertex and projection factor (—G)”(|)”, and the color factor 
are taken into account at the end of the calculation. We can relate iLn+i to Kn by one-gluon 
crossing, the effect of which can be expressed by an operator TZ: 


Kn+l=TZ[Kn]. 


(32) 


Now let us expand Kn as 


16 

Kn — Oin ' b — ^ ^ Oin,lbj , 
I=l 


(33) 


where an,i are coefficients of each basis bj. For example, Ki = — 7 ^ (g 7 ^, so that 017 = —fes/. 
The operation of TZ on each basis 6/ is expressed as 

P-ibi] = - (7m'S'+^^/ g lp.S-li) , (34) 

where the factor (—1) comes from the charge conjugation explained in Eq. (17). Dehning the 
recursion matrix R by"^^ 


Rjibj = TZ[bi\, (36) 

the recursion relation between Kn and Kn+i gives 

Oin+1,1 = RijOin,j- (37) 

Note that the matrix R is a function of 7l(s/4), B(s/4) and s = P^. The explicit form of the 
matrix Rjj is shown in Table 2. The coefficient vector becomes 

Otn = R^~^CXi, (38) 

In the practical calculation, we define an inner product I among bj’s as 

I(bi, bj) = T:rD[uiUjsliljs\ + TrD[uiUjsljlis\ (35) 

with s = 1 + 27 ®, to invert Eq. (36) to obtain Rij. Note that the choice of s is arbitrary as long as the rank of 
the matrix I{bi,bj) is 16. 
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Figure 6: Diagrammatic expression of the resummation of the completely crossed ladder diagrams 
performed in subsection 3.3. 



Figure 7: Values of the determinant Eq. (11) obtained by the resummation of the completely crossed 
ladder diagrams performed in subsection 3.3 for various current quark masses. Each line corresponds 
to the current quark mass toq = 0.012, 0.1, 0.2, 0.5, 1.0 x rj [GeV]. 


and thus the sum K of the completely crossed diagrams is given by 

OO OO 

K = Y,Kn = (i)” 

n=\ 
oo 

n=\ 

= g\( [(-G) (I) i? + 6] - [(-G) (f) 4i? - 3] 


n=l 


Ql 


b, 


(39) 


where we used the color factor shown in Eq. (44). 


4 Numerical results and discussion 

After the resummation of the completely crossed ladder diagrams, the diquark BS eq. actually 
has a bound-state solution. Figure 7 is a plot of det(/ — H{P‘^)) in the scalar diquark channel 
as a function of for various current quark masses mo. The figure shows that the BS eq. 
does have solutions of det(7 — H) = 0. This is the main result of our study. 

We observe similar behavior in the axial-vector diquark channel, and the resulting masses 
are summarized in Figure 8. The axial vector diquarks tend to be heavier than the scalar 
diquarks, which is consistent with other studies [23,24,44-46]. 

There are two subtleties in our calculation that have to be discussed. One is the ambiguity 
of the truncation scheme for the quark self-energy, as we have mentioned repeatedly in the 
previous section. We check that the existence of bound-state solutions is unchanged under 
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Table 2: Recursion matrix Rij. 


1 2 3 4 5 6 7 8 



1 

2 

iAB 

-A?sl2 

-2B^ 

—iABs 






3 

A 

-iAB 

-2B^ 

-A^sl2 

iABs 






5 

-B^ 

-iABs/2 

iABsl2 

-Ah^A 






6 

iAB 12 

B^ 

J^sjA 

-iABsl2 

iAB 

-^2 s/2 

-2^2 

—iABs 


7 

-iABjl 

J^sjA 

B^ 

iABs j 2 

-iAB 

-2^2 

-A^sl2 

iABs 

T 

8 

-A?IA 

iAB 12 

-iABl2 

-B^ 






9 





-B^ 

-iABsl2 

iABs/2 

-A^s'^/A 


10 





iAB 12 

^2 

A^s/A 

-iABs/2 


11 





-iABl2 

^2s/4 

S2 

iABs/2 


12 





-J^jA 

iAB 12 

-iABl2 

-B^ 


13 










14 










15 










16 




J 





Rij 








9 

10 

11 

12 

13 

14 

15 

16 


1 

2 

3 

4 

5 





4^12^ 

—8iABs 

8iABs 

16S2s 


6 





8iAB 

—AA^s 

-16^2 

—8iABs 


7 





-8iAB 

-16^2 

— 4^2s 

8iABs 


8 





im'^js 

8iAB 

-8iAB 

AA?s 

1 

9 

-4^2 

-2iABs 

2iABs 







10 

iAB 

-A?sl2 

-2B^ 

—iABs 

2iAB 

4^2 

A?s 

-2iABs 


11 

-iAB 

-2B^ 

-A?sl2 

iABs 

-2iAB 

^2s 

4^2 

2iABs 


12 





-A^ 

2iAB 

-2iAB 

-4^2 


13 

-B^ 

-iABs/2 

iABsl2 

-Ah^A 

—A?s 

2iABs 

-2iABs 

-AB^s 


14 

iAB 12 

B^ 

J^sjA 

-iABsl2 

-iAB 

A?sl2 

2^2 

iABs 


15 

-iABl2 

A^sjA 

B^ 

iABsl2 

iAB 

2^2 

A^s/2 

—iABs 


16 

-A?IA 

iAB 12 

-iABl2 

-B^ 

-AB‘^/s 

-2iAB 

2iAB 

—A'^s 


12 
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Figure 8: Diquark masses as a function of the input current quark mass. These values are obtained 
by the resummation of the completely crossed ladder diagrams performed in subsection 3.3. 

different truncation schemes for the DS equation, although the values of bound-state masses 
change. The results in Figure 7 and Figure 8 are calculated with the quark self-energy Eq. (6), 
since it is the leading term in the ’t Hooft limit. The other is the space-like poles appearing 
when the current quark mass is uiq < 0.1 x r/ [GeV]. This can be a consequence of the 
simplification of the QCD model and the truncation scheme we adopt, and there is a possibility 
that other models or truncation schemes resolve this artifact. 

We would like to note that the presence of the diquark bound state does not contradict 
its confined nature. In the previous studies [19,21,22] the absence of diquark bound states is 
regarded as a realization of confinement. This argument is based on a fact that correlation 
functions without one-particle poles do not contribute to the S-matrix. However, there is 
another realization of confinement based on the BRST quartet mechanism [47], in which it 
is proven that under a certain condition color non-singlet states are also BRST non-singlet 
states, and hence they do not appear in the physical Hilbert space. Therefore, if quarks are 
confined by this mechanism, then diquarks are confined simultaneously. 

Our result that the diquark BS eq. has a bound-state solution gives a natural explanation 
to phenomenologically successful diquark constituent models: first two of three (or four) quarks 
form a bound state and can be treated as a one-particle state. However, we should note that 
the absence of bound-state solutions to the diquark BS eq. itself does not contradict to these 
models. What we have investigated in this paper (and in other papers about the diquark BS 
eq.) is merely a pure two quark system, after all. There is a possibility that more complex 
dynamics happening inside hadrons makes two quarks bound even if their pure BS eq. has no 
bound-state solutions. 

We also comment on our interpretation of the completely crossed ladders. In weakly 
bound systems like the hydrogen atom, the constituent particles (proton and electron) can 
be regarded as nearly on-shell, and the leading effect of the exchange of photon becomes 
instantaneous interaction. In this case, the completely crossed diagrams of higher orders 
require the constituent particles to become off-shell by the energy-momentum conservation, 
which leads to suppression of these effects [48]. However, in the case of diquarks and mesons, 
a naive application of above discussion is impossible because the on-shellness of quarks inside 
them are no longer meaningful. This interpretation is valid if the current masses of quarks are 
small. The reason why the coupling expansion of the meson BS kernel is successful may be 
the suppression of the color factor discussed around Eqs. (22)-(25). To summarize, we suppose 
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that the completely crossed ladders can be neglected when the constituents are regarded as 
nearly on-shell (hydrogen atom) or when their coefficients are suppressed (meson), and diquark 
is an exception of both cases. 


5 Conclusion 

The diquark Bethe-Salpeter equation with the Munczek-Nemirovsky model is investigated. We 
find that the contributions of completely crossed ladder diagrams in the diquark channel affect 
the behavior of the Bethe-Salpeter kernel significantly, in contrast to the meson channel. We 
perform a resummation of all the completely crossed ladder diagrams in the Bethe-Salpeter 
kernel. As a result, we find that diquark bound-state solutions exist. Although the values 
of resulting diquark masses may contain model- and truncation-artifacts and are far from 
a reliable prediction at this stage, the presence of the bound-state solution is rather stable 
against variations in the parameter values and truncation schemes. 

Acknowledgements 

The work of R.J. and G.M. was supported by Grant-in-Aid for JSPS Fellows (No. 25-8360 and 
No. 26-10887, respectively). The authors are grateful to the referee for the careful reading of 
the manuscript and a number of excellent suggestions. 

A Color Factors 

The basic formula is 


{t°')ij{t°‘)ki = \ {5ii5kj — jfSijSki) 


(40) 


where (t“)ij is the generator of SU(A^) in fundamental representation. In the meson channel, 
the ladder diagrams have a factor 



2 " \ ' N) — 1 


(41) 


whereas the completely crossed ladder diagrams have 


^ ^ ^ ^ (43) 


The ladder diagrams in the diquark channel have a factor 


^^^2 , , , -j-Ojn—l 


-TanjTan-i _ _ _ ^Ta2^Tai 



(43) 


whereas the completely crossed ladder diagrams have 
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Table 3: Coefficients and 


a( 

n) 

j 


1 

2 

3 

4 

5 

6 


1 

2 







2 

-8 

-12 





n 

3 

32 

64 

16 




4 

-128 

-336 

-256 

-112 




5 

512 

1536 

1728 

1152 

192 



6 

-2048 

-7360 

-12288 

-12672 

-6144 

-1472 



n) 

j 


1 

2 

3 

4 

5 

6 


1 

4 







2 

-12 

-8 





n 

3 

32 

32 

16 




4 

-112 

-256 

-336 

-128 




5 

384 

1152 

1728 

768 

64 



6 

-1472 

-6144 

-12672 

-12288 

-7360 

-2048 


B Practical expressions of the DS and BS eqs. in the coupling 
expansion 

Here we present concrete expression of equations used in Section 3. The DS eq. with the 
quark self-energy (26) becomes a pair of self-consistency equations, 

A ^ 

= 1 - (t)" 4“' E (45) 

n j=l 

VD ^ 

- E(-o)” E (46) 

n j=l 

where s = with ^(s) and H(s), and the coefficients are summarized in Table 3. 

The solutions of these equations are plotted in Figure 1 for the cases of up to 0(G), O(G^) 
and 0(0^). 

In Section 3, we see that the problem of whether a homogenous BS eq. has nontrivial 
solutions or not reduces to a problem of whether a determinant equation (11) has solutions or 
not. The matrices we used in that section are expressed as 

4f(») = E(-G)”(i)"<!!.“>fc!.“'(«) (47) 


in the meson case, and 


H{s) = Y,i-Gr {lTd^n^h(^\s) (48) 

n 
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in the diquark case. Here 


and 


,(M) 


,(M) 


u{M) _ 'V(l.l) 

- \ ^(M) ^(M) I - 


n,(2,l) n,(2,2) 


,(D) 


.(D) 


= I ftCtf)'"’ I = E 


h'- 

n,(2,l) n,(2,2) 


.K 


h 

(M)' 


(M) 


fc=0 \ n,(2,l)j^ 


h 


h 

(D) 


(D) 

72,(1,1),/C 


k = 0 \ 72,(2,1),/c 


h 


(M) 
n,(l,2),fc 


{B/A) 


{A/B) 


h 




h 


(D) 

n,(l,2),fc 


{A/B) 


h 


(D) 

n,{2,2),k 


{Ah)'^-’^{B‘^)’‘, (49) 


(M) 

n,{2,2),k / 


{BIA)\ 


(50) 




where s = A = ^(s/4) and B = i3(s/4). The coefficients and are 

summarized in Table 4. The values of det(/ — H) calculated using the above matrices are 
plotted in Figure 4. 
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Table 4 : Coefficients in Eqs . ( 49 ) - ( 50 ). 


^(M) 




k 





k 

0 

1 

2 

3 

4 

5 

6 

0 

1 

2 

3 

4 

5 

6 

n 

(bj) 








n 

(hJ) 









( 1 , 1 ) 

1 

-4 







( 1 , 1 ) 

-1 

4 






1 

( 1 , 2 ) 

4 

0 






1 

( 1 , 2 ) 

-4 

0 







( 2 , 1 ) 

0 

2 







( 2 , 1 ) 

0 

-2 







( 2 , 2 ) 

- 1/2 

2 







( 2 , 2 ) 

1/2 

-2 







( 1 , 1 ) 

- 1/4 

3 

8 






( 1 , 1 ) 

- 1/4 

3 

8 





9 

( 1 , 2 ) 

-2 

-4 

0 





9 

( 1 , 2 ) 

-2 

-4 

0 






( 2 , 1 ) 

0 

-1 

-8 






( 2 , 1 ) 

0 

-1 

-8 






( 2 , 2 ) 

1/2 

3 

-4 






( 2 , 2 ) 

1/2 

3 

-4 






( 1 , 1 ) 

1 

6 

12 

-16 





( 1 , 1 ) 

- 1/4 

0 

0 

16 




o 

( 1 , 2 ) 

3 

20 

48 

0 




Q 

( 1 , 2 ) 

- 3/2 

-8 

-24 

0 





( 2 , 1 ) 

0 

3 

28 

48 





( 2 , 1 ) 

0 

- 3/2 

-16 

-24 





( 2 , 2 ) 

- 1/2 

0 

36 

80 





( 2 , 2 ) 

1/2 

3 

-12 

-32 





( 1 , 1 ) 

- 13/16 

- 23/4 

-19 

28 

128 




( 1 , 1 ) 

- 1/16 

1/4 

5 

76 

128 



A 

( 1 , 2 ) 

- 5/2 

-23 

-112 

-176 

0 



A. 

( 1 , 2 ) 

-1 

-11 

-64 

-80 

0 




( 2 , 1 ) 

0 

- 11/4 

-28 

-92 

-160 




( 2 , 1 ) 

0 

- 5/4 

-16 

-44 

-64 




( 2 , 2 ) 

1/2 

7/4 

-19 

-92 

-208 




( 2 , 2 ) 

1/2 

19/4 

5 

4 

-16 




( 1 , 1 ) 

1 

33/2 

99 

204 

144 

-64 



( 1 , 1 ) 

- 1/16 

- 3/2 

-9 

36 

96 

64 



( 1 , 2 ) 

11/4 

45 

288 

720 

704 

0 



( 1 , 2 ) 

- 7/8 

-15 

-108 

-240 

-224 

0 



( 2 , 1 ) 

0 

13/4 

51 

288 

816 

832 



( 2 , 1 ) 

0 

- 11/8 

-21 

-108 

-336 

-352 



( 2 , 2 ) 

- 1/2 

- 3/2 

51 

396 

1296 

1472 



( 2 , 2 ) 

1/2 

21/4 

9 

-36 

-336 

-512 



( 1 , 1 ) 

- 61/64 

- 159/16 

- 93/2 

-78 

-36 

720 

2048 


( 1 , 1 ) 

- 1/64 

- 27/16 

- 21/2 

66 

492 

1680 

2048 


( 1 , 2 ) 

- 21/8 

- 139/4 

-240 

-936 

-2528 

-3264 

0 


( 1 , 2 ) 

- 3/4 

- 73/4 

-168 

-648 

-1472 

-1344 

0 


( 2 , 1 ) 

0 

- 51/16 

- 79/2 

-234 

-960 

-2224 

-2688 


( 2 , 1 ) 

0 

- 21/16 

-23 

-162 

-672 

-1168 

-768 


( 2 , 2 ) 

1/2 

45/16 

- 9/4 

-78 

-744 

-2544 

-3904 


( 2 , 2 ) 

1/2 

105/16 

123/4 

66 

-168 

-432 

-64 


























C Comment on the number of the independent elements of 
the basis 

The following four elements made of ^ and 7 ’s in Eq. (31) 

liJilvlplu ® Iplulpla , Iplvlpla ® flplulpla , 

flplvlplu ® Iplvlpla, flplvlplu ® tlplvlplu, 

may seem independent of the other 16 elements in Eq. (31), but these are not since the following 
equation holds: 

Iplulpla «) Iplvlpla = - ® tip + ^Iplv ® Iplv + -^tlplulp ® tlplulp- (51) 

To check this, note that Eq. (51) is rewritten as 

- 4P“P^/p2^ ilalplulp < 8 ) Ifilplulp - ^lalp < 8 ) l/3lp) = 0. (52) 

This is equivalent to 

_ 4p«p/3/p2^ {lalplulpM-ffi-i^'y^-ip - AjalpMji^jp) = 0. (53) 

for an arbitrary Hermitian 4x4 matrix M, which can be expanded by ipi---lp,„- One can 
readily check that the second parentheses of Eq. (53) vanishes when n is odd, while it is 
proportional to 5ap when n is even and vanishes after it is multiplied by the first parenthesis. 
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